home *** CD-ROM | disk | FTP | other *** search
/ Graphics Plus / Graphics Plus.iso / general / fractal / kaos.lha / fixptlib / fode_user.c < prev    next >
Encoding:
C/C++ Source or Header  |  1989-11-18  |  622 b   |  35 lines

  1. /*
  2. ### compute a periodic orbit of period r for vector fields ###
  3.  
  4. Input:    r: 0: equilibrium point, >0: periodic orbits
  5.     time:    time
  6.     time_step: (r=0) time_step for a next orbit
  7.     x[n]:    x(time)
  8. Output: x[n]:    r=0: value of a vetor field, f(x)
  9.         r>0: r'th image of an orbit x(time+period*r)    
  10. */
  11.  
  12. fode_user(x,time,r,dim)
  13. int r,dim;
  14. double x[],time;
  15. {
  16.     int i,j;
  17.     extern int model;
  18.     extern double *t_v,*param;
  19.     extern int (*f_p)();
  20.  
  21.  
  22.     if(r==0){
  23.         (int) f_p(t_v,0,x,param,time,dim);
  24.         for(j=0;j<dim;j++) x[j] = t_v[j];
  25.     }
  26.     else {
  27.         /*
  28.         for(i=0;i<r;i++){
  29.             fode_one_user(t_v,x);
  30.             for(j=0;j<dim;j++) x[j] = t_v[j];
  31.         }
  32.         */
  33.     }
  34. }
  35.